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A Bayesian analysis is carried out to identify the consistent regions of the mSUGRA parameter 
space, where the newly- discovered Higgs boson's mass is used as a constraint, along with other 
experimental constraints. It is found that mi/ 2 can lie in the sub-TeV region, Aq/tuq is mostly 
confined to a narrow strip with |Ao/rao| < 1, while mo is typically a TeV or larger. Further, 
the Bayesian analysis is used to set 95% CL lower bounds on sparticle masses. Additionally, it 
is shown that the spin independent neutralino-proton cross section lies just beyond the reach of 
the current sensitivity but within the projected sensitivity of the SuperCDMS-lT and XENON-IT 
experiments, which explains why dark matter has thus far not been detected. The light sparticle 
spectrum relevant for the discovery of supersymmetry at the LHC are seen to be the gluino, the 
chargino and the stop with the gluino and the chargino as the most likely candidates. 
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Introduction: The most recent search at the LHC [TJ [2] 
for the Higgs boson [3 with the combined 7 TeV and 
8 TeV data indicates a signal for the Higgs boson with 
mass 125.3 ± 0.4 (stat.) ± 0.5 (syst.) GeV for CMS with 
a local significance of 5.0<j and with mass 126.0 ± 
0.4 (stat.) ± 0.4 (syst.) GeV for ATLAS with local signifi- 
cance of 5.9 a. As is well known the Higgs boson mass at 
the tree level lies below the Z° boson mass, but it can be 
made larger by inclusion of loop corrections. However, 
in supergravity grand unification [4 there is another up- 
per limit, i.e., of about 130 GeV due to the constraint 
of radiative breaking of the electroweak symmetry (for a 
review see [5 ) as well as other experimental constraints 
(for a recent analysis see [6] [7]). The correction to the 
Higgs boson mass is given by [5] (for a review see [9]) 
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where X t = A t — \i cot /3, where A t is the Aq parameter 
run down to the weak scale (see Eq. ([2])), v = 246 GeV, 
and Ms is an average stop mass. The loop correction in 
Eq. |l]) is maximized when X t ~ \^6Ms. There are also 
additional loop corrections from, e.g., the 6-quark sector 
as well as from higher loops. The early searches at the 
LHC-7 gave a possible hint of the Higgs boson in the mass 
range ~ (117— 129) GeV [10] and the combined Tevatron 
analysis reported an excess between (115 — 140) GeV [TT] . 
These findings have led to significant activity [6] [I2j [13] 
to investigate the implications of the results for super- 
symmetry. 

Implications for mSUGRA: We note that the scale Ms 
in Eq. ([!]) which is determined by the soft parameters 
depends sensitively on the Higgs mass. In the analysis we 
use the Higgs boson mass constraint within the Bayesian 
statistical framework to estimate the soft parameters of 
mSUGRA (sometimes referred to as CMSSM) which are 
given by [4 



where mo is the universal scalar mass, m^ is the uni- 
versal gaugino mass, Aq is the trilinear couplings and 
tan/3 is the ratio of the two Higgs VEVs in MSSM, 
and fi is the Higgs mixing parameter. The soft param- 
eters of Eq. define our model's parameter set, = 
{mo, mi/2, Ao, tan and additionally we consider a set 
of the most sensitive standard model nuisance param- 
eters, i/j = ^m t ,m b (m b ) MS ,a s (m z ) MS ,a EM (m z ) Ms y 

These together form the basis parameter set: = 
Using Bayes's theorem, the posterior probability density 
function (PDF) for the theory described by O, which may 
be mapped to observables, £(0) to be compared against 
experimental data, d is given by: 
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where C = p(d\^(Q)) is the likelihood function-the terms 
of which are described in Table [ij 7r(0) is the distri- 
bution in prior to considering experimental results, 
and Z = p(d) is the Bayesian evidence which can be 
used in model selection. However, in our goal of param- 
eter estimation, it serves only as a normalization fac- 
tor. We present results obtained by considering both 
the 2D marginalized posterior PDF (where the full N- 
dimensional posterior PDF of Eq. ([3| has been integrated 
over the other parameters), as well as the profile like- 
lihoods (where the confidence levels are determined by 
comparison to the global best-fit point). (For a more 
detailed description see the second reference in [15].) 

The analysis was done by using SusyKit [19], which 
employs the MultiNest [18] package for sampling pa- 
rameter points efficiently, and uses SoftSUSY [20] for 
spectrum calculation, and micrOMEGAs [21] to cal- 
culate the relic density as well as for the indirect con- 
straints. The credible intervals, marginalized posterior 
PDF's, and profile likelihood distributions were calcu- 
lated using the plotting routines of SuperBayes [22] . 
which is largely based on the tools provided by Cos- 
moMC [23]. The constraint from the — 2 measure- 
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TABLE I: Summary of the observables used to estimate the 
mSUGRA parameters. Only the upper-half of the Gaussian is 
used in the consideration of Q/i 2 , i.e., there is only a penalty 
for values larger than the central value which allows for mul- 
ticomponent dark matter |30j . The 95% CL limits have been 
evaluated under the assumption of only theoretical uncer- 
tainty, so the distribution used here is based on the error 
function, given explicitly in the fourth reference of [22] . 



merit is not imposed in this analysis and this issue will 
be discussed later in the text. 

In our analysis, we took our prior knowledge of the 
parameters to be either flat linear distributions or flat 
logarithmic distributions, with mo G (0.05, 8) TeV (log), 
mi/a G (0.05, 5) TeV (log), A G (-30, 30) TeV (linear), 
and tan/3 G (3,60) (linear). We have fixed sgn(/i) to be 
positive. The Standard Model nuisance parameters were 
allowed to vary in 2 a windows of their central values, as 
quoted in Table [I] Our MultiNest sampling parameters, 
as defined in [18 , were nii ve = 20,000 and tol=0.0001. 
It has been shown in [24] and the second reference in [18] 
that these parameters are not only sufficient to provide 
a map of the posterior PDF, but also to find the true 
best-fit point which is essential for the profile likelihood 
analysis. 

In our likelihood analysis we use the CMS result since 
that result was available earlier pQ. We report our fits 
to the data, including the Higgs mass, in Fig. [I] in the 
form of 2D posterior PDF maps (left panels) as well as 
the profile likelihood maps (right panels). The posterior 
mean is marked with a large dot and the global best-fit 
is marked with a circled 'X'. (Note that while the best- 
fit point is crucial in Frequentist likelihood-ratio tests, it 
has no significance in the Bayesian framework.) The top 
panels exhibit the constraint in the mo — mi/2 plane and 
show that mo is typically a TeV or larger, while mi/ 2 can 



lie below 500 GeV. The middle panels exhibit the con- 
straint in the Aq/tuq — tan/3 plane, and here one finds 
that most of the allowed parameter space lies in the nar- 
row strip l^o/mol < 1 with a small strip in the range 
l^o/mol G (—2,-6). The bottom panels exhibit the con- 
straint in the tua — tan f3 plane, and here one finds that 
the majority of the allowed range of tua lies above 1 TeV. 
Thus rriA ^ rn h o for the majority of the parameter space 
and thus we are in the so-called decoupling limit. 

It was pointed out in [25] that most of the experimen- 
tally consistent parameter space of mSUGRA lies on the 
Hyperbolic Branch (HB) [26 , 27 of radiative breaking of 
the electroweak symmetry under the LHC-7 constraints. 
The HB region has sub-regions which we may label as Fo- 
cal Point (HB/FP), Focal Curves (HB/FCi, i=l,2), and 
Focal Surfaces (HB/FS). It was shown in [6] [25] that the 
HB/FP is mostly depleted while the remaining parame- 
ter space lies on HB/FCi or HB/FS. Specifically we note 
here that the right edge of Aq/ttiq in Fig. [I] is ~ 1. The 
value \Aq /mo | = 1 was argued as string- motivated in [28] 
and was shown to be the asymptotic limit on the focal 
curve HB/FCI in [25]. 

In Fig. [2] we present the 2D posterior PDF's (left pan- 
els) and the profile likelihoods (right panels) in the planes 
of the phenomenologically important sparticle masses. 
The top panels present the results in the gluino-squark 
mass plane, and indicate that the gluino can be below a 
TeV. The second row is plotted in the squark-chargino 
mass plane and demonstrates that the chargino masses 
are only bounded from below by the direct searches at 
LEP. The next row exhibits our fit in the stau-stop mass 
plane. Here one finds that the stau and stop masses are 
typically large except for a small strip where the stop 
mass can lie below a TeV. This is largely to be expected 
as we rely on a heavy stop to provide a sizable loop cor- 
rection to the Higgs mass. The bottom panels show the 
analysis in the fi — rrig plane. One finds that fi is typically 
quite light, i.e., ji can be significantly below 500 GeV. 

Using the marginalized ID posterior PDF we are able 
to set lower limits on the sparticle masses from the 2<j 
credible regions. We present those limits here: nig > 
1.39 TeV, m~± > 196 GeV, rriA ~ ttih ~ mH± > 
1.3 TeV, m h > 3.1 TeV, m fl > 3.1 TeV, m~ q > 5 TeV, 
and m-g > 4.8 TeV. The profile likelihood analysis yields 
different results. Here, we find the 95% CL sparti- 

690 GeV, m~± > 95 GeV, 
ir > 580 GeV, m?, > 



cle lower limits to be m~ g > 

mA ~H ~H± > 540 GeV, : 



310 GeV, rriq > 1.5 TeV, and m l > 580 GeV. We note 
that as expected the lower limits given by the profile like- 
lihood analysis lie lower than the limits given by the PDF 
analysis. The analysis thus indicates that the light parti- 
cles in mSUGRA in view of the Higgs mass measurement 
are the neutralino, the chargino, the gluino, the stau and 
the stop. Among these the most likely candidates for 
discovery in the next phase of CERN experiment are the 
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FIG. 1: Left panels: plots of the 2D posterior probability 
densities, 1 a and 2 a contours are also drawn. Right panels: 
plots of the profile likelihoods. Top: in the mo — mi/2 plane. 
Middle: in the Aq/ttiq — tan/3 plane. Bottom: in the vtla — 
tan/3 plane. The posterior mean is marked by a large dot 
while the best-fit point is shown by a circled 'X'. The color 
bar above the top panel gives the relative likelihood which 
increases left-to-right. 



gluino, the chargino and the stop. 

125 GeV Higgs boson and dark matter: Neutralino- 
proton spin independent cross section cr?o depends sen- 

XiP 

sitively on the Higgs boson mass (for a discussion see 
[6]). Thus considering the ~ 125 GeV Higgs mass leads 
to a more constrained prediction for dark matter. In 
Fig. |3j we give a plot of 1Z x cr?o p as a function of the 
lightest neutralino mass m^o where the factor 1Z is de- 

IS 



fined by K = (Qh 2 ) / (Vth 2 



'WMAP< 



and (yih 2 



'WMAP 



the central value of the WMAP-7 data. By only apply- 
ing a likelihood penalty for points that are above the 
WMAP-7 limit, we have taken into account the possi- 
bility that there may be additional components of dark 
matter beyond the neutralino [30] . Quite remarkably, the 
bulk of the credible region of mSUGRA falls essentially 
exclusively between the current limits on dark matter by 
XENON-100 [31. and the projected sensitivity of Super- 
CDMS [33] and XENON-IT [32]. 

We discuss now the constraint from — 2. In su- 
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FIG. 2: Left panels: plots of the 2D posterior probability 
densities, 1 a and 2 a contours are also drawn. Right panels: 



plots of the profile likelihoods. Top: in the nig 



plane. 
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l x plane. Bottom: in the u- 



Lower-middle: in the 
g plane. The posterior 



mean is marked by a large dot while the best-fit point is shown 
by a circled 'X'. The color bar above the top panel gives the 
relative likelihood which increases left-to-right. 



persymmetric theories, sparticle loops make significant 
contributions to the anomalous magnetic moment of 
the muon [29] if the relevant sparticles (charginos, 
neutralinos, smuons, sneutrinos) entering the loops 
are relatively light. The experimental determination 
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FIG. 3: (color online) Plots of K = (Qh 2 ) / (ft/i 2 ) WMAp vs 
the neutralino mass rn^o. The left panel presents the 2D 
posterior PDF, and the right panel presents the profile like- 
lihood. The analysis shows that virtually all of credible re- 
gion of mSUGRA will be probed by the SuperCDMS [33] and 
XENON- IT [32] experiments. The color bar above the panels 
gives the relative likelihood which increases left-to-right. 



pends sensitively on the hadronic correction to the 
standard model value. There are two main proce- 
dures for the estimation of the hadronic correction, 
which are either using the e + e~ annihilation cross 
section or from r decay. The result using the e + e~ 
annihilation gives Sa^ = (28.7 ± 8.0) x 1(T 10 (3.6 a) 
while for r-based hadronic contributions one has 
Sa^ = (19.5 ± 8.3) x 1(T 10 (2.4 a) [34 . In any case, 
within the universal soft SUSY-breaking paradigm there 
would be tension between the — 2 result (specifically 
the one using e + e~ annihilation cross section) and 
the 125 GeV Higgs boson mass since the mo scale is 
rather high. If the — 2 results stay, then there are at 
least two avenues open to make compatible the g^ — 2 
results and the Higgs boson mass. The first possibility 
is that we stay within the universal soft breaking 
paradigm and additional contributions to the Higgs 
mass arise due to the presence of extra matter which 
can generate new loop corrections to the Higgs mass, 
or from extra gauge groups under which the Higgs is 
charged yielding corrections to the Higgs mass through 
D-terms. Alternatively, one could give up universality of 
soft parameters and consider non-universal or flavored 
SUGRA models [35]. For instance, to satisfy the g^ — 2 
constraint one may consider the soft scalar mass for 
the first two generations much smaller than for the 
third generation, or the sleptons being lighter than the 
squarks. These possibilities require further investigation. 

Conclusion: In this work we have analyzed the impli- 
cations of the Higgs boson discovery at CERN for super- 
symmetry. Specifically we analyzed the mSUGRA model 
to delineate constraints on soft parameters and identified 
the light particles that are prime candidates for discov- 
ery in the next phase of runs at the LHC. The analysis 
presented here explains why supersymmetric dark matter 
has not been seen thus far since essentially all of the pa- 



rameter space lies below the current sensitivity of dark 
matter experiments due to the high Higgs mass. The 
analysis also points to excellent prospects for the discov- 
ery of dark matter at SuperCDMS and XENON- IT as 
well as the possibility of light neutralinos, charginos and 
gluinos, and possibly light stops and staus. 
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